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Abstract 

At future linear e''~e^ collider experiments in the TeV range, Sudakov dou- 
^ ■ ble logarithms originating from massive boson exchange can lead to significant 

corrections to the cross sections of the observable processes. These effects are 
important for the high precision objectives of the Next Linear Collider. We use 
the infrared evolution equation, based on a gauge invariant dispersive method, 
to obtain double logarithmic asymptotics of scattering amplitudes and discuss 
how it can be applied, in the case of broken gauge symmetry, to the Standard 
Model of electroweak processes. We discuss the double logarithmic effects to 
both non-radiative processes and to processes accompanied by soft gauge boson 
emission. In all cases the Sudakov double logarithms are found to exponentiate. 
We also discuss double logarithmic effects of a non-Sudakov type which appear 
in Regge-like processes. 
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1 Introduction 

The Next Linear Collider (NLC) will explore e~^e~ processes in the TeV energy regime, 
and probe the Standard Model of elementary particles to great accuracy. Electroweak 
processes which may reveal New Physics, such as supersymmetry, are especially inter- 
esting. Therefore the accurate calculation of the scattering amplitudes of such high 
energy processes in the Standard Model is very important. The main corrections to 
the Born amplitudes at high energies, ^/s, are double logarithmic (DL) contributions 
of the form (g'^log (s/m?))^ which arise from soft vector boson exchanges. Here we 
derive the DL asymptotics of scattering amplitudes, for a sequence of gauge theories 
leading up to the Standard electroweak Model. 

The DL correction for QED processes dates back many years to the paper by 
Sudakov [^. In Section 2 we present the original Sudakov form factor for QED in a 
form in which it can be generalized to non-Abelian gauge theories, a subject to which 
we then turn. We start with the processes where all the scalar products pipj of the 
momenta of participating particles are of the same order. In Section 3 we compute the 
DL corrections to processes governed by massless non-Abelian gauge theories, and then, 
in Section 4, we turn to the interesting case of broken gauge symmetries, discussing, 
in particular, the DL effects for Standard Model electroweak processes. 

In all these cases the DL terms can be resummed in exponential forms which are 
generalizations of 

M = MBorneXp Ll f^ iy,(s, /i^)") . (1) 

Here Wi is the probability of the emission of a soft and quasi-coUinear gauge boson 
from an external particle i, and the summation is over all n external charged particles 
participating in the process. The virtual particles are subject to an infrared cut-off fi 
on their transverse momenta. 

We use an equation for the evolution of the scattering amplitudes as a function of 
the infrared cut-off /x of the transverse momenta of virtual particles. Our approach is 
based on a gauge invariant dispersive method which is a generalization of a powerful 
theorem on photon Bremsstrahlung due to Gribov 0. His remarkable theorem, proved 
by dispersive methods, states that, due to gauge invariance, the region of applicabil- 
ity of well known formulas for accompanying photon Bremsstrahlung is considerably 
extended at high energies. In these formulas the amplitude for a process with the emis- 
sion of a soft photon is expressed, in a factorized form, in terms of the amplitude for 
the non-radiative process, with the particles taken on-mass-shell. Since the amplitude 
is taken on-shell, gauge invariance is guaranteed. The evolution equation approach 
greatly simplifies the computation of DL effects for electroweak processes, which are 
mediated by broken symmetry with the photon having components in both the SU{2) 
and U{1) gauge groups. 

The evolution equations for the amplitude as a function of the infrared cut-off 
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parameter /z are analogous to the Renormalization Group Equations (RGE). Basi- 
cally we start from the domain of very large /i, where the Born amplitude applies, 
and evolve down to small fi, matching the solution at the mass thresholds bounding 
every new kinematic domain. It is therefore not surprising that the structure of the 
exponentiation of the DL effects has a clear physical interpretation. 

Of course for observable processes it is necessary to consider the DL corrections to 
both the non-radiative and the radiative processes. For radiative processes we again 
have exponentiation of the Sudakov DL corrections in a form similar to (|l|). The 
dependence on the infrared cut-off parameter n is canceled in the measurable semi- 
inclusive process, and is replaced by a dependence on the experimental acceptance cuts 
on the soft boson emissions. In Sections 3.1, 4.1 and 4.2 we address these issues. 

As well as DL corrections of the Sudakov type, there are DL contributions specific 
to forward or backward scattering, when a final particle is quasi-collinear to an initial 
particle and their energies are almost equal. This is the domain of Regge kinemat- 
ics and is the subject of Section 5. Here we have the DL corrections related to the 
exchange of soft pairs of fermions or of bosons in the crossed channel. The infrared 
evolution equations are now non-linear in the amplitude^, and in the j-plane repre- 
sentation are of Ricatti-type form. In simple cases the latter forms can be reduced 
to a Schrodinger-type equation with a harmonic potential. This approach can also 
be applied to production processes in the multi-Regge kinematics. In this case the 
evolution equations are solved in a sequence of domains starting from the region where 
the Born amplitude applies and evolving down in the cut-off /x to the region where soft 
particles are emitted. The vector boson and fermion reggeization can be also easily 
verified in the DL approximation. Moreover, the infrared evolution equation allows 
the construction of scattering amphtudes with quasi-elastic unitarity. 

2 Sudakov form factor in QED 

The high energy asymptotics of electromagnetic processes was calculated many years 
ago within the framework of QED [jl[-|p. In particular the amplitude for e^e~ elastic 
scattering at a fixed angle [s ~ \t\ ~ \u\ ^ m? 3> A^, where m is the electron and A a 
fictitious photon mass) in the DL approximation has the form 



M=MBornr'[^,^] , (2) 



.2 I J ^ 

m?'' A^ 



^An analogy is now to the non-linear RGE for the coupling constant, whereas for the 
Bremsstrahlung processes an analogy is the linear RGE for the moments of structure functions. 

^A plays the role of the infrared cut-off. In physical cross sections the divergence in A of the elastic 
amplitude is canceled with the analogous divergences in processes with soft photon emissions. 



where A^Bom is the Born amphtude for e^e scattering and V is the Sudakov form 
factor. The DL approximation apphes in the energy regime 
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s s rn 

— ~ a log ^ log — 



alog ^ ~alog^log-— ~ 1, (3) 



where the QED coupling a = c^/Att <^ 1. Thus each charged external particle effec- 
tively contributes vT to the total amplitude. The Sudakov form factor appears in the 
elastic scattering of an electron off an external field |l|]. It is of the form: 





To specify R it is convenient to use the Sudakov parametrization of the momentum of 
the exchanged virtual photon : 

k = vpi + up2 + k± , (5) 

where pi and p2 are the initial and final momenta of the scattered electron. R{s) can 
then be written as the integral over u and v. 

' du f dv f ^ \ , ) f ^ \ , \ e{suv - X') , (6) 

where s ~ |t| ~ 2piP2- The first two factors in the integrand correspond to the propa- 
gators of the virtual fermions which occur in the one-loop triangle Sudakov diagram. 
The 9 - function appears as a result of the integration of the propagator of the photon 
over its transverse momentum fc^, noting that the main contribution comes from the 
region near the photon mass shell [||]: 

SUV = \^ + k\. (7) 

To DL accuracy (6) gives for \ <^m: 

R 

where the result comes equally from two different kinematical regions, v ^ u and 
M 3> f . Therefore one can write R = 2r. 

We can obtain physical insight by presenting the two equal contributions separately. 
In the first region, with v ^ u, the virtual photon is emitted along pi and the parameter 
V is given by the ratio of energies of the photon and the initial electron. Here instead 
of u, it is convenient to use Eq. (|^) to replace it by the square of the transverse 
momentum component of the photon. Then integrating over v and fcj^ gives 
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dv r"'"^ dki f dki /•min(|fcx|/m, 1) (j^y 



A/v^ V JA2 kj_+m'^v'^ J\^ fc]_ J|fcx|/v^ V 



(9) 



in the DL approximation, which may be evaluated to give half of R. The quantity r is 
proportional to the probability Wi of the emission of a soft and almost collinear photon 
from an external particle with energy ^/s and mass rrii, i.e. 



w,{s,y) = -r[—,-^. (10) 

If several charged particles participate in a process, for example e^e' — > f fff , then 
analogous contributions appear for each external line, provided all external invariants 
are large and of the same order. This leads to the general result 

A^ = -MBornexp(--^w;i(s,A2)j , (11) 

where n is the number of external lines corresponding to charged particles. In sum- 
mary the soft emissions described by the Sudakov form factor is a quasi-classical effect 
which does not depend on the hard dynamics of the process. In particular there are 
no quantum mechanical interference effects in the DL Sudakov corrections, for large 
scattering angles. 

3 Generalization to non-Abelian gauge theories 

Sudakov effects have been widely discussed for non-Abelian gauge theories, such as 
SU{N) and can be calculated in various ways (see, for instance, 0). We consider 
here the scattering amplitude in the simplest kinematics when all its invariants sij = 
2pipj are large and of the same order sij ~ s. A general method of finding the DL 
asymptotics (not only of the Sudakov type) is based on the infrared evolution equations 
describing the dependence of the amplitudes on the infrared cutoff fi of the virtual 
particle transverse momenta [0]. This cutoff plays the same role as A in QED, but, 
unlike A, it is not necessary that it vanishes and it may take an arbitrary value. It can 
be introduced in a gauge invariant way by working, for instance, in a finite phase space 
volume in the transverse direction with linear size / ~ 1/yU. Instead of calculating 
asymptotics of particular Feynman diagrams and summing these asymptotics for a 
process with n external lines it is convenient to extract the virtual particle with the 
smallest value of \k±\ in such a way, that the transverse momenta |fc'_i_| of the other 
virtual particles are much bigger 

k'l > fci > /i^ . (12) 

For the other particles fc^ plays the role of the initial infrared cut-off /i^. 

In particular, the Sudakov DL corrections are related to the exchange of soft gauge 
bosons, see Fig. 1. For this case the integral over the momentum k of the soft (i.e. 




Figure 1: Feynman diagrams contributing to the infrared evolution equation ( p!3|) for 
a process with n external legs. In a general covariant gauge the virtual gluon with 
the smallest value of k± is attached to different external lines. The inner scattering 
amphtude is assumed to be on the mass shell. 

|/i;°| <^ \/s) virtual boson with the smallest fc_L can be factored off, which leads to the 
following infrared evolution equation: 

2 (27r)4 j,i^j^i -''*»fci»A'2 k^ + «e [kpj){kpi) 
xT\j)T\l)M{pr,...,Pn.kl), (13) 

where the amplitude A^(pi, ...,p„; fc_|_) on the right hand side is to be taken on the 
mass shell, but with the substituted infrared cutoff: /i^ — > k\. The generator 
T"'{T){a = 1, ...yN) acts on the color indices of the particle with momentum pi. The 
non-Abelian gauge coupling is g. In Eq. (p!3D, and below, k± denotes the component 
of the gauge boson momentum A; transverse to the particle emitting this boson. Note 
that in Sudakov DL corrections there are no interference effects, so that we can talk 
about the emission (and absorption) of a gauge boson by a definite (external) particle, 
namely by a particle with momentum almost coUinear to k. It can be expressed in 
invariant form as A;_|_ = iam((kpi){kpj) / (pipj)) for all j ^ I. The above factorization 
is related to a non-Abelian generalization of the Gribov theorem^ for the amplitude of 
the Bremsstrahlung of a photon with small transverse momentum k± in high energy 
hadron scattering p. 

The form in which we present Eq. (|T3|) corresponds to a covariant gauge for 
the gluon with momentum k. Formally this expression can be written in a gauge 

■^The non-Abelian generalization of Gribov's theorem is given in ( pl| ) below, together with a de- 
scription of its essential content. 



invariant way if we include in the sum the term with j = I (which does not give a 
DL contribution). Indeed, in this case we can substitute PiPj by —p^p^^d^y^k), where 
the polarization matrices of the boson d^jj{k) in the various gauges differ by the terms 
proportional to k^ or k'^ giving a vanishing contribution due to the conservation of the 
total color charge Y^aT"" = 0. Thus we have the possibility of choosing appropriate 
gauges for each kinematical region of quasi-coUinearity of k and pi. We can, however, 
use (|13|) as well, noting that in this region for j y^ I we have PjPi/kpj ~ Ei/uj, where 
El is the energy of the particle with momentum pi and uj the frequency of the emitted 
gauge boson, so that: 



Mipi,...,Pn;iJ'^) = MBorniPl, ■■■,Pn)- -7^;^ J2 J 



2^2 n ^^^ .min(|fcx|M, 1) dy 



X CiM{px,...,Pn\k\) , (14) 

where Ci is the eigenvalue of the Casimir operator T°-{l)T°-{l) {Ci = Ca for gauge bosons 
in the adjoint representation of the gauge group SU{N) and Ci = Cp for fermions in 
the fundamental representation). 

The differential form of the infrared evolution equation follows immediately from 



dM{pu...,Pn;fJ'^) 



aiog(/i2) 

where 



K{f/)M{pi,...,Pn;fJ^'^), (15) 



with 



W.isy)^£,C.r[^,f). (17) 

As in the Abelian case, Wi is the probability to emit a soft and almost coUinear 
gauge boson from the particle / with mass mi, subject to the infrared cut-off /i on the 
transverse momentum. Note again that the cut-off fi is not taken to zero. The function 
r is determined by (9) for arbitrary values of the ratio rrii/^. To logarithmic accuracy, 
we obtain from (0): 
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The infrared evolution equation (|T5|) should be solved with an appropriate initial con- 
dition. In the case of large scattering angles, if we choose the cut-off to be the large 
scale s then clearly there are no Sudakov corrections. The initial condition is therefore 

M{pi, ...,Pn; S) = MBorniPl, ■■■,Pn), (19) 



and the solution of (|T5[) is thus given by the product of the Born amphtude and the 
Sudakov form factors: 

A<(pi,...,p„;/i^) = A<Born(Pl,...,Pn)eXp f -- X^ W/(s, /i^) j (20) 

Therefore we obtain an exactly analogous Sudakov exponentiation for the gauge group 
SU{N) to that for the Abelian case, see ([l^). Theories with semi-simple gauge groups 
can be considered in a similar way. 

3.1 DL corrections to processes with soft emissions 

Since ultimately we are interested in measurable cross sections we have to consider the 
DL corrections to amplitudes of processes with soft emissions, as well as those without. 
Only in inclusive cross sections will the dependence on the infrared cut-off /i^ disappear, 
being replaced by parameters specifying the experimental acceptance. To put it in 
another way, cross sections of the emission processes receive large (DL) contributions 
from regions where the emitted bosons are soft and the emission angles are small. 
Therefore, to be consistent, we need to calculate cross sections of such processes as 
well. This is easy to do for QED processes, where the single gauge boson (the photon) 
is neutral and does not possess self-interactions. Therefore soft photons are emitted 
independently according to a Poisson distribution. In non-Abelian theories, gauge 
bosons are not neutral and interact with each other. Consequently, the soft emission 
does not follow a Poisson distribution p[. 

We again consider the simplest situation, when the additional soft gauge boson is 
emitted in the process with all invariants sij large. Of course, for the emission of a 
boson almost collinear to the particle the direction of the particle with momentum pi, 
the invariant 2kpi is small in comparison with s. In the case of non-Abelian gauge 
theories the corresponding amplitude for the emission of a soft gauge boson with small 
k^ <^ jj? has, according to the Gribov theorem, the following form: 

j=l '^Pj 

The possible corrections to this factorized expression are of the order of k^/ jj?. How- 
ever, to DL accuracy, we can substitute /x^ in the arguments of the scattering am- 
plitudes by its boundary value k\. Notice that the amplitude on the r.h.s. of ( pT]) 
is taken on-the-mass shell, which guarantees its gauge invariance. The result (^) is 
highly non-trivial in the Feynman diagram approach. It means, that the region of ap- 
plicability of the classical formulas for the Bremsstrahlung amplitudes is significantly 
enlarged at high energies. V. Gribov proved this result by using dispersion relations in 
the variables {k -|-pj)^, and demonstrating that for small |fc_|_| the pole terms in these 





Figure 2: Schematic Feynman diagrams contributing to the real gauge boson emission 
for a process with n external legs. For fcj^ <C fi^ the diagram on the left corresponds 
to the non-Abelian generalization of Gribov's theorem (|2T|). The diagram on the right 
leads to an additional term, (p2D, in the kernel of the evolution equation in the case 
when k^ :» /x^. 



invariants are much larger than the corresponding cut contributions due to the gauge 
invariance of the theory p[. The region of applicability of (|21|) corresponds to the 
situation when the momentum of the emitted soft boson does not spoil the kinematics 
of the non-radiative process. This implies that the frequency uj of the boson emitted 
off the particle with momentum pi should be much smaller than the energy Ei, and 
that the emission angle ■(?; ~ |fc_L|/ci; should be much smaller than the angle between 
pi and all other momenta pj (otherwise one needs to include interference effects). For 
the physical Coulomb gauge {eq = 0, e ■ k = 0) in the kinematical region where the 
gauge boson is emitted at a small angle "i^j with respect to the particle with momentum 
Pi in ( pll) only the term with j = i contributes. The method based on the infrared 
evolution equation allows us to calculate, in the DL approximation, the amplitudes of 
the hard processes accompanied by the emission of any number of soft gauge bosons 
0. Let us consider the amplitude for the emission of one soft gauge boson in the 
region of quasi-coUinearity of its momentum with momentum pi (i.e., the emission of 
a soft boson by a particle with momentum pi). When the transverse momentum |fc_|_| 
of this boson is much less than the infrared cut-off /i, used for the virtual particles, the 
amplitude is given by the term with j = i in the sum in (^). But we need to know the 
emission amplitude in the opposite case, \k±\ :» /i. It can be found from the evolution 
equation in this region using expression (PH) at /i = \kj_\ as the initial condition. The 
kernel of the evolution equation in this region differs from the corresponding kernel in 
the region \k_j_\ -C /i (that is the kernel (16) of the evolution equation for the amplitude 
without emission) by a term connected with the emission of a virtual boson from the 



real gauge boson with momentum k: 






where I^a is given by ( [LTD with C; = C^i and m^ = //, see Fig. 2. It is clear, that 
this new term in the kernel for evolution from fcj^ to /i^ leads to an additional term 
WA{kj_, fi^) in the Sudakov exponential. Thus, the amplitude for the emission of one 
gauge boson with small transverse momentum from the hard scattering process is of 
the form: 

j=l '^Pj 

/ 1 " 1 \ 

xexpi--J2Wi{s,fi')--WA{kl,fi')j . (23) 

We note again that here fc^ means the square of the component of the three-dimensional 
momentum transverse to the momentum of the emitting particle, say for example, pi. 
We can write fc^ in the invariant form (kpi){kpj)/{pipj) with j ^ I, which does not 
depend on j. 

Again we see that we have the exponentiation of the Sudakov DL corrections. 
Note that in the Abelian QED case we have Wa = and the exponent for the photon 
Bremsstrahlung amplitude remains the same as that for the process without photon 
emission. It is related to the Poisson distribution for soft photon production. 

The exponentiation of virtual DL corrections holds for multiple emission processes 
as well. In QED it is trivial, since the soft photons are emitted independently. In non- 
Abelian gauge theories it is not so simple. The main complexity is connected with the 
nontrivial structure of the amplitudes for multiple emission of real soft gauge bosons, 
arising from their self-interaction. But in the DL approximation these amplitudes can 
be calculated. Due to the coherence effect, the branching cascade develops only in the 
region of sequentially shrinking angular cones |^. In this region the Born amplitudes 
for multiple emission processes have a factorized form and the virtual corrections expo- 
nentiate [§] . It is proved by solving the infrared evolution equation in a series of regions 
where the infrared cut-off /x is bounded between a sequence of decreasing transverse 
momenta of the emitted gluons, using in each region the solution of the previous region 
as the initial condition. The final result is 

Mipi, ...,p„; fcl, ..., K] ^J?) = A<Born(Pl, ■■■,Pn; h, ••■, K) 
/ 1 " 1 ^ \ 

^^p[-^j:wi{s,f.')--j:wA{ki,t,')j , (24) 

where fcj are the momenta of the emitted gluons with strongly ordered energies and 
ki± are their components transverse to the momenta of the emitting jets. 
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4 Sudakov effects in broken gauge theories 

The same method, based on the infrared evolution equation, is also applicable to broken 
gauge theories. Let us consider for definiteness the Standard electroweak theory, where 
the physical gauge bosons are a massless photon (described by the field A^) and massive 
W^ and Z bosons (described correspondingly by fields Wj^ and Z^). To DL accuracy, 
all masses can be set equal: 

M^ ~ M,4/ ~ Mniggs ~ M 

and the energy considered to be much larger, ^/s ^ M. In the unbroken phase 
the corresponding Abelian and Yang-Mills fields are denoted by B^ and W^^, with 
a = 1,2, 3. The physical fields are linear combinations of the fields in the unbroken 
theory with coefficients depending on the Weinberg angle 6^. The left and right handed 
fermions are correspondingly doublets (T = 1/2) and singlets (T = 0) of the SU{2) 
weak isospin group and have hypercharge Y related to the electric charge Q, measured 
in units of the proton charge, by the Gell-Mann-Nishijima formula Q = T^ + Y/2. 

In the evolution equation in the DL approximation the value of the infrared cutoff 
fi can be chosen in two different ranges : 1) ^/s ^ ^^ M and 2) /i -C M. The second 
case is universal in the sense that it does not depend on details of the electroweak 
theory. It will be discussed below. In the first region we can neglect spontaneous 
symmetry breaking effects, in particular gauge boson masses, and consider the evo- 
lution equation in the unbroken phase with effectively massless particles B and W"". 
Of course one could calculate everything also in terms of the physical fields A'^ ^ Z^ 
and W^ . In the unbroken phase this is equivalent to the description in terms of the 
original fields B^ and W^ and leads to the same final result, but the intermediate steps 
will be more complicated because there are cancellations between non-exponentiating 
terms from diagrams with Z and 7 exchanges. The separation of these contributions is 
not gauge invariant and if we would consider the diagrams without virtual photons we 
would violate SU{2) x f/(l) symmetry. Taking into account only such diagrams leads 
to nonexponentiating DL effects in an axial gauge §]. Fig. 3 illustrates this at the 
two loop level. The loss of gauge invariance is related to the fact that the photon field 
contains the component Wj^ of the non- Abelian field W^}, and so omitting the virtual 
photons would violate the conservation of the weak isospin current (in the unbroken 
theory). In region 1) the infrared evolution equation, written in terms of the unbroken 
fields, is of a form analogous to (|15D if we assume, for simplicity, that all the charged 
particles have masses rrti < M 



dM{pi,...,Pn;iJ.) log(s//i 



91og(/i2) (47r) 



^ /T,(T, + 1)+^' 




Here Tj is the total weak isospin of particle i, Yi is its weak hypercharge, and g and 
g' = gtanO^ are the couplings of the SU{2) and U{1) gauge groups, respectively. The 
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Figure 3: Two- loop 'rainbow' Feynman diagrams contributing to DL corrections in an 
axial gauge. The photon contribution has DL corrections in both the regions fi^ <^ M^ 
and M^ ^ yU^ ^ s. Taken together with the W- and Z-contributions, it yields the 
exponentiation of the Sudakov DL terms in the electroweak theory. In the region 
M^ <^ /i^, the spontaneously broken gauge symmetry is restored and omitting the 
photon contributions would lead to a non-gauge invariant result. 
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sum in (^) is to be performed over all n external particles. As before, the initial 
condition is given by the requirement that for the infrared cut-off /i^ = s we obtain 
the Born amplitude. The solution of ( ]25|) is thus given by 



M{pi 



1 Pn; A'- 



M 
X exp 



Born(Pl, •••,Pn) 

2^„ /,.2 



i^|:H.,.,./(| 



i=l 



(26) 



The expression in the brackets in the exponential can be written in terms of the pa- 
rameters of the broken theory as follows: 



j'i 



Yi 



^ = e' + ^' (^^(^^ + 1) - (^')') + ^ir (^' - ^i^' ^-^0 



V 2 / ^ ^ cos^ du, 

where the three terms on the r.h.s. correspond to the contributions of the soft photon 
(interacting with the electric charge Cj = Qig'sva.Q^^ the W^ and the Z bosons, re- 
spectively. Although we may rewrite solution (^) in terms of the parameters of the 
broken theory in the form of a product of three exponents corresponding to the ex- 
changes of photons, W^ and Z bosons, it would be wrong to identify the contributions 
of the diagrams without virtual photons with this expression for the particular case 
el = 0. This becomes evident when we note that if we were to omit photon lines then 
the result would depend on the choice of gauge, and therefore be unphysical. Only for 
9u, = 0, where the photon coincides with the B gauge boson, would the identification 
of the ef term with the contribution of the diagrams with photons be correct. 

Choosing the cutoff /x in the second region, /i <^ M, the infrared evolution equation 
takes the following form: 



dM{pi,...,Pn;iJ'^) 



E 



e? 



log 



max(m?,/x2) 



M{pi,...,pn,l^^). (27) 



51og(/i2) f^, (47r)2 

Evidently, only the photon contribution remains in this region. Now the appropriate 
initial condition is given by (|2^) evaluated at the matching point ^ = M. The solution 
is thus 



M{pi 



,Pn;fJ' 



MBovniP 

X exp 



Bornl/^l, ■■■,Pn) 



2(47r) 



X: r,(T, + l)+tan2^^M 



i=l 



^'^'w^, 



1 n 



X exp 

^ 4 = 1 

A^Born(Pl,---,Pn) 
-.2 



X exp 



X exp 
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2(47r) 

n / 

■El 

i=l 



Y.(\Ti{T, + l)+tan''e^(^ 



i=l 
^2 



^°^'ii^, 



{Any 



log 



M' 



m,M 



log ^ + log ^ log 



mt 



mf 



m,- 



/^^ 



,(2^ 
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where the last equahty holds for fi <^ rrii and mf <^ M^ from the respective expansions 
of Wi in ([ToD . Let us stress that (^) and (|28|) are applicable for processes involving 
chiral fermions as well as gauge bosons, provided that all the invariants are large {0{s)) 
compared to M^. Note, that in the case, when quarks or gluons participate in the 
reaction, we should multiply these expressions by the Sudakov factors corresponding to 
the virtual gluons emitted by these colored particles. The infrared evolution equations 



([25| ) and (|27|) have a clear physical meaning analogous to the renormalization group 
equations and therefore it is natural to expect that the next-to-leading corrections to 
the kernels can be calculated. 

For physical observables soft real photon emission must be taken into account in an 
inclusive way and effectively the parameter /i^ in (pSD will be replaced by parameters 
depending on the experimental requirements. 

4.1 DL effects for electroweak processes with soft emission 

The calculation of amplitudes for processes with the emission of a gauge boson in the 
kinematical region which gives DL contributions to the cross sections (i.e. the region of 
soft quasi-collinear emission) is similar to the analogous calculation for unbroken gauge 
theories, with complications of the type that we discussed above. One is that we have 
to consider separately two regions of |fc_L| of the emitted boson: first \k±\ <^ M and 
second |fc_L| ^ M. At high energies the cross sections of the emission processes receive 
DL contributions from both of these regions. Of course, W^ and Z boson emissions 
contribute in the second region only. Therefore consideration of the first case is very 
simple. For values of the infrared cut-off jj? ^ k^ the amplitude for the emission of a 
soft photon by a particle with momenta pi (i.e. emission within the cone along pi not 
containing the momenta of the other particles) has, in the physical (Coulomb) gauge, 
a factorized form 

S*Vi 

M{pi,...,pn;k;ix^) = ei - — M{pi,...,Pn;iJ^^), (29) 

kpi 



according to the Gribov theorem of (21). However, if kj^ <^ M^, the kernel of the 
infrared evolution equation does not change when the cut-off /i^ changes from the 
domain /x^ ^ kj^ to fj? -C fcj^ (since the W^ and Z bosons do not contribute in this 
first region). Therefore, ( pUj ) remains valid at arbitrary values of the cut-off fi. 

In the second region, |fc_|_| ^ M, the result is more involved. We need to start 
from fi^ ^ k\, where we can use a generalization of Gribov's theorem. Consider again 
the case when a gauge boson is emitted by a particle with momentum pi. We have 

M«(pi, ...,p„; k; /x^) = GS(0|^ M{pu ...,p„; fx') , (30) 
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where M.{pi, ...,p„;/i^) is given by (^^ and 



G^ = —;=T^ for W^ emission, 

° V2 

G^ = — ^ (t^ - Q sin^ 9^) for Z emission, 
cost^^ ^ ^ 

GQ=Qgsiia9w for 7 emission, (31) 

with Q = {T^ + V/2). Then we have to solve the evolution equation in the region 
M^ <^ IJ? <^k^ with the initial condition given by (|30|) at the matching point jj? = k^. 
In fact, it is more appropriate to work in terms of the fields W^ and By. The kernel of 
the evolution equation remains unchanged for the emission of the S-particle due to its 
Abelian nature. On the other hand the emission of the W"" -particle leads to the same 
additional contribution as in the unbroken theory (see (p2D ). Therefore in the cut-off 
region fc^ ^ /i^ ^ M^, we obtain 

^«(pi, ...,p„; k- /i^) = G?(/)^ A^(pi, ...,p„; /i^) , (32) 

kpi 

with the same amplitude A^(pi, ■■■,Pn', A^^) as in (p^) and 

Of = G^ + gco^e^T^(e^v[-\wA{kl,^?)^-l^ , 

Gl = Gl + gs^ne^T^[e^^[-UvA{kl,^^'')^-lY (33) 

Finally we study the region of the infrared cut-off /i^ <^ M^. In this region the kernel of 
the evolution equation is determined by the electromagnetic interaction only; therefore, 
the only contribution related to the emitted particles is that for H^^emission. This 
contribution is given by 

^^'^(^ ) = "2 91og(/.^) = (i^^^^ii^ ('') 

where ww{k\, ^'^) is defined by (|1^) and (9) with M^ = M^. Consequently, for values 
of the infrared cut-off /x^ <C M^ we obtain: 

>('^(pi,...,p„;fc;/i2) = G«(/)^A<(pi,...,p„;/i2), (35) 

kpi 

where A^(pi, ...,Pn; /i^) is given by (^) and 

G± = GU^^(^-]^WA{kl,M^)-]^ww{kl,i^) + Uvw{kl,M^)Y 
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G^ = Gl + g sin 9^T^(^exp(^-^WA{kl,M^)^-iy (36) 

The important point here is the difference of the DL exponent for a nonradiative 
process and the process with photon emission, which leads to a violation of the Poisson 
distribution for photons in the DL approximation at high energies. It is a direct 
consequence of the fact that the photon has a non-Abehan component. 

4.2 DL effects in semi-inclusive cross sections 

Measurable cross sections have an inclusive nature (at least with respect to photons, 
since only cross sections with an infinite number of emitted soft photons are observ- 
able). Let us consider the DL corrections to such cross sections. 

It is clear that for the emission of real gauge bosons the same cut-off /x^ must 
be used as for virtual ones. Therefore, to calculate an experimentally measured cross 
section we have to take the cut-off /i^ less than the lower bound nl^p of k^ of those 
photons emitted in processes which are not included in the cross section. 

The calculation of the cross section is simple if the experimental conditions are such 
that only processes with emission of photons with fc_|_ < M^ are allowed. In this case 
the non-Abelian component of the photon is not essential, so that photon emissions 
obey a Poisson distribution. Therefore for the cross section with an arbitrary number of 
emitted photons with momenta lying inside regions Qi of the momentum space around 
the emitting particles with momenta Pi, we obtain: 

d(r{pi, . . . ,p„) = do-elastic (Pl, ■■■,Pn) exp^W^^p) , 

where rfcXeiastic is the cross section of the non-radiative process and w^^p is the probabihty 
of the emission of photons with fc^ > /i^ inside the allowed region 






W^exp 



n „2 ^ j2 



= ^hl i:^rr^4^i7^-^(fcl-/^')- (37) 

^ Atx'^ Jn, kj_ + m^ujyEf to 

Since the upper bound on k^ of the photons which are allowed to be radiated is less 
than M^, we must use the cut-off /i^ < M^ and, consequently, ( [28| ) for the matrix 
element of the non-radiative process. Therefore, we obtain 



da{pi,...,Pn) = daBorn{Pl,---,Pn) 

X exp 



^' tiUm + l)+tanHJ^)']\og' ' 



15 



xexp -J2Q'^[w,{s,fi^)-w,{s,M^))+w2,^ , (38) 

where Wi{s, fi"^) is given by (10) and (9). Evidently, the dependence on /x in X) Qfwi^s^ /i^) 
and ly^xp cancels in the exponential. 

If we include in the observed cross section the emission of gauge bosons with 
transverse momenta larger than M, then the problem becomes much more complicated 
because of the non-Poisson distribution of soft emission in non-Abelian gauge theories 
0. Let us consider here the simplest example of the cross section completely inclusive 
of photons emission in the two-loop approximation. Then the cross section can be 
written as 

da{pi, ...,Pn) = dcTBoTniPl, ■■■,Pn) \l + Sv + ^+Sr + Srv + ^\ , (39) 

where 6y is the one-loop virtual correction, 6r comes from one-photon emission taken in 
the Born approximation and 6rv from the one-loop correction to one-photon emission. 
Due to exponentiation of the DL terms of the Sudakov-type in virtual corrections, 
the term 6^/2 in ( |39D gives the two- loop virtual correction. The term S^/2 gives 
the correction from two-photon emission in the Born approximation, since in this 
approximation the two photons are emitted independently. 

In the considered case of the cross section completely inclusive of photon emission 
the cut-off parameter /i^ can be taken as large as M^ (but not greater, because the 
cross section does not include W^ and Z emission). Each of the corrections considered 
above depends on yU^, but their sum in (|39|) does not; therefore, we can take the most 
suitable value of the cut-off to calculate the cross section. It is easy to see that the 
most convenient choice is /i^ = M^. In this case from (26) we have: 

^.=-ilji|:({^.ffi+i)+'-^».(|)]>o^^]^) ■ («) 

The correction due to one-photon emission taken in the Born approximation is 6r = 
w2xp where w^^p is given by (|37| ) with /x^ = M^, Ei ~ y/s and the region Qi defined by 
inequality u < Ei. It gives 

S. = ± Q?^.(., M^) =tQlj^. log^ ^ ■ (41) 

The one-loop contribution to this correction Srv can be found with help of (26), (31), 
(32) and (33): 

f^Air^JM'^ kj_ J\k^\ uj \- 
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E 



e2 



■^.g^iog^^-g,7;3^iog^^ 



(42) 



- (47r)^ 

The important point here is that the virtual correction to the one-photon emission does 
not coincide with the corresponding correction to the non-radiative process (which 
means violation of the Poisson distribution) and depends on the momentum of the 
emitted photon. It is a consequence of the fact that the photon has a non-Abelian 
component. 

One can check that using an arbitrary cut-off /i^ < M^ gives the same result, but 
the calculation is more complicated. 

5 DL effects for amplitudes with Regge kinematics 

For two particle scattering in either the forward s ^ (— t) or the backward s ^ 
(— m) directions, corresponding to the Regge regime, there are situations, where the 
amplitudes have DL contributions different from those of the Sudakov type. These 
logarithms appear from Feynman diagrams in which the soft particles with the minimal 
transverse momenta are pairs of virtual fermions or bosons exchanged in the crossed 
t or M channels. For QED and QCD the infrared evolution equations for such Regge 
processes are known [4, 5, 7]. The new features presented by the Standard Model of 
weak interactions are that the gauge group is semi-simple, SU{2) x t/(l), and that there 
is a large difference in the particle mass scales. The general strategy of finding DL 
asymptotics in this case is to first solve the evolution equation in the region /i ^ Mi, 
where Mi is the largest particle mass in the theory; then to solve it in the region 
M2 -C /i <^ Ml, where M2 is the next largest particle mass, using the solution of the 
previous equation at /i = Mi as the initial condition, and to proceed with these steps 
until the intermediate particles are light quarks, electrons, neutrinos and photons. 

For simplicity we consider only the evolution equation in the region where the cut- 
off parameter // is much larger than all particle masses. In this case it is natural to 
calculate the amplitudes in terms of the massless particles of the unbroken theory, that 
is the leptons and quarks, the Higgs, B and W"" bosons. By solving the generalized 
infrared evolution equation in the DL approximation we sum up contributions of ladder- 
type diagrams in the crossed channel, with all possible insertions of the gauge bosons, 
leading to the Sudakov double logarithms. Because the transverse momenta of the 
virtual particles are strongly ordered, the integral kernels in the ladder diagrams are 



given by the splitting kernels of the DGLAP evolution equ ation [|Ty] which describes 
all possible transitions among fermions and bosons in the Standard Model. These 
splitting kernels can be simplified since the Regge kinematic domain corresponds to a 
strong ordering of the Bjorken variables Xi along the ladder. 

For definiteness we study the simple case of the backward lepton (/) - antilepton 
(/) scattering for s ^ — t ^ — m. In the Born approximation, the contribution of the 
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Feynman diagrams with the Higgs, photon, B and W"" boson exchanges depend on 
the hehcities Q = Xi/2 (with Aj = ±1) of the initial and final fermions, which have 
momenta 

Pi, Ph Pi' - Ph Pv - Pi- (43) 

We consider the scattering of the leptons e and z/g belonging to the first generation 
with the smallest masses. Because the coupling of the Higgs particle to leptons is 
proportional to their mass, we can neglect the Higgs contribution. The Dirac spinors 
describing the / and I states with the definite helicities A are 

".(rt=/w(t;J. «.(-P) = \/w (!-,;_ J . (44) 

where 

-— -V5a = Av9a- (45) 

\P\ 

With the use of these expressions we can calculate the matrix elements of the 7-matrix 

structures which appear in the Born approximation for backward lepton-antilepton 

scattering, ll ^1 1', 

-vxj{-pjhaUxXPi)ux^,iPi'hcTVx-i-PiT) = as' ' '' '' =25a,,-a/v,-a-<5a,,-v(46) 

-y^A^(-Pr)7at^A^(-PF)^v(PF)7aMA,(pO = a/ ' '' '' =26x^,XjSx,,,x-Sxi,x,,- (47) 

It is important, to note that both the helicity structures of a^ and at are proportional 
to Sx, A- , which means, that the helicities of the two fermions in the crossed u channel 
are opposite in sign. Therefore the Born diagram with a virtual W"' boson gives a 
negligible contribution to the backward scattering. Because the B boson interacts 
with hypercharge Y, its contribution to the Born amplitude is also zero for the right- 
handed neutrinos (if they were to exist). 

Below we consider only the non-trivial case, when only e^ , i^- and !/_,_ can partic- 
ipate in the reaction. Because the hypercharge Y = —1 for left-handed fermions and 
Y = —2 for the right-handed electron, in this case the Born amplitude is given by 

m'^^T" '" = y (a'' '"'" '" + a'' '''" '") . (48) 

In the DL approximation we also obtain a significant simplification of the helicity 
structure of the scattering amplitude: 

M'^'^'^''- = ^(a^''^'"'-/(.//.^)+ar^'"'"/(-^//i^)), (49) 

fis/fx') = /+(s//i2) + r(s//x2), (50) 
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where the functions /^(s//i^) describe contributions with positive and negative signa- 
ture: 

f^i-s/fi') = ±/^(3//.^), 

fLnis/f^') = 1, /Born(^//^') = 0. (51) 

The amphtudes /^ are assumed to satisfy dispersion relations of the form 

fHs/f^') = -- r (^ ± ^) Im/^(.7^^)^ , (52) 

71 Jo \S — S S + S' / S 

where in the Born approximation liaf^^^^{x) is different from zero only in the region 
where x is of the order of unity. In particular, from the known behavior of the Born 
contribution /Born(s//i^) at large s and its degeneracy in the signature we obtain 



--rim/B^_(.V/^^)^ 
vr Jo s 



The absence of degeneracy of the amplitude with respect to signature is a result of the 
presence of the non-planar diagrams which arise from the virtual soft 5-bosons in the 
infrared evolution equations (cf. the QCD case 0): 



rflog(/i2) 

log+(-V/"') 


g'2 .log 

Stt^ Jo 
(3g' + 9g" 


87r2 

dlogifx^) 
log+(-^^') 


V 4 

g'2 .log 

STT^io 

f3g' + 9g" 



r^ fHs,/^?) f r'=Mm/+(.2//i^)rflog^) rflog^+ 
Jo \Jo U'^ u^ 



\og-i-s/fi') fSg' + g'^' 
87r2 



fHs/f^') + ^^^^^-4^ --^L_LJL_ /-(,/^2)^ (53) 



/"'""/-(V/^^) /°'"lm/-(V/.^)rflog^ rflog^+ 
Jo \Jo w^ w^ 



f-l I 2.. log (-g//i^) ('ig^+jf^ ,+ , I 2n .P-.N 

where we introduced the notation 

log^(-s/;.^) = ^-^^-'1^"^^ l"g(^//^') . 

The linear terms in j^^ on the right-hand side of the evolution equations, are 
the Sudakov contributions which arise from soft W"" and B virtual bosons coupling 
to external lines with different momenta pi c^ pjj and pu — Pj- The non-linear terms 
appear due to soft pairs of leptons exchanged in the crossed u channel. 

Note that, in the above equations, the term related to the s-channel Sudakov vertex 
with W^ exchange is proportional to the Casimir operator T(T+ 1) = 3/4 for the weak 
isospin group SU{2), and the total B boson contribution from the four Sudakov vertices 
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having singularities in the s and t channels is proportional to (—1 — 2)^/4 = 9/4. The 
corresponding contribution of these vertices to the terms mixing the amplitudes with 
different signatures is proportional to 3/4 for W^ exchange and (—1 + 2)^/4 = 1/4 
for B exchange. Non-linear terms take into account the contribution of the ladder 
diagrams for two leptons interacting through B exchanges. 

The above equations can be simplified if we use the Mellin representation for the 
amplitudes: 

where C^{j) are the signature factors 

exp(— iTTJ) ± 1 



eu) 



sin(7rj) 



and where the integration contour L is situated along the imaginary axes to the right 



± 



of all singularities of the t-channel partial waves /,• . In the j-representation, taking 



into account the degeneracy of the partial waves with different signatures in the Born 
approximation, the evolution equations can be written in the form of Ricatti equations 

[7]: 

■'^ 167r2j2V^iy ^^ 327r2 J jdj \ 327i^ J f' ^ ^ 

if we take into account only the terms with effective DL variables 

1 g'^ 1 2g'^ 



Lo^ IGvr^ j2 ' 3(j(2 -|_ g(j('2 



(58) 



and neglect other terms. Furthermore, using the following definition for the t-channel 
partial waves 

V^^M = /f, (59) 

the equations can be written as follows 



LU I au duj a uP- 



V? (CJ) = 1 - ( ^-^ \ - -- ^ + ^— . (61) 
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Both equations can be reduced to linear form after introducing new functions ip{uj) 
and x(^) according to the definitions 

(/9+(c<j) dip{(jj)/duj uj V^~(^) dx(}jj)/duj oj 1 

uj a%p{uj) 2 ' CO ttxi^) 2 2auj 

They have the form of Schrodinger equations 



d 



2 ^2 



^^, + -.^ + a^ + -j^M=0, (63) 



V rfcu2 ' 4 4^2 ' ' V ' V^ V ^(^) 

The first is an equation for the harmonic oscillator. Therefore taking into account the 
boundary condition 

lim ■i/'fcj) ~ exp ( —10"^ I cu"' , (65) 

ui^+oo \4 / 

so as to match onto perturbation theory, we obtain the solution 

^{cu) = DaiV^co) , (66) 

where Dp{z) is the parabolic cylinder function: 

D,{z) = '.^^tm r e—^x'^-^dx . (67) 

T{-p) Jo 

Therefore for the scattering amplitude with positive signature we have, to DL accuracy, 

f*''^"''^ = Lil i^J i|log(exp(-«.V4)fl.(V^.)). (68) 

with UJ = Anj/g'. The equation for negative signatures can be easily solved for a = 
—1/4 (corresponding to complex values of the Weinberg angle): 

X{uj) = v^ exp(-cuVl6) ^ Q, 1, y j , (69) 

where \E'(q;, 7, z) is the confluent hypergeometric function: 

^(a, 7, z) = —— / e-""a;"-Vl + xV-^'^^dx . (70) 

T[a) Jo 

Therefore for a = —1/4, the amplitude with negative signature, in the DL approxima- 
tion, is 

ris/^') = Zllf^ (l.y i, d log ^ fl, 1, ^] , (71) 

■' ^ '^ ' 2 Jl 27: \fiy a-' d] ^ \8' ' 8 J ' ^ ' 
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with uo = AiTJ/g'. For general values of the parameter a the solution can be obtained as 
a perturbation series in uj~^. This perturbative expansion is convergent for all values 
of the DL parameter g"^ \og^ {s / fi^) . 

Because the function Dp{z) does not have any singularities in the complex plane, 
the high energy asymptotics of the scattering amplitude with positive signature 

r(^//i^)~3^ J^ = J-^ (72) 

47r \/—a 



is determined by the position of its zeroes, Dp[zi) = 0, in the left-half plane, Rezi < 0. 
The function "^{a, 1, z) has a singularity ~ Inz as z — > 0, and therefore the scattering 
amplitude with negative signature behaves as 

ris/fi')^\n-'s (73) 

at high energies. Such behavior of /^(s//i^) is valid also for other values of the param- 
eter a. Note, that similar DL asymptotics were obtained in QED for e^e~ backward 
scattering amplitudes. For energies smaller than the masses of W"" and Z bosons, only 
diagrams with virtual photons give rise to a DL contribution. The photons interact 
with electric charge Q = T^ + Y/2, which is the same for left- and right-handed elec- 
trons. It is interesting, that at very large energies, the behavior of e~^e~ backward 
scattering amplitudes is different from the QED prediction even in the case of the 
Abelian Standard Model with g = 0. 

The DL asymptotics of lepton-antilepton forward scattering amplitudes, which is 
related to t-channel exchange of two fermions and an arbitrary number of W"" and B 
bosons, can be obtained in a similar way. In this case one should take into account the 
diagram with a virtual W"' boson even in the Born approximation. The reggeization 
of the fermions and W"' boson can also be verified in the DL approximation with the 
use of the infrared evolution equation. Moreover, one can construct amplitudes with 
quasi-elastic unitarity |ll|, which is important for the theory of the BFKL Pomeron 



| T2| . Some other applications of the DL asymptotics in QCD are given in [|T3[, |T^ . 



6 Discussion 

We have calculated the double logarithmic (DL) corrections to the amplitudes of a 
number of different high energy processes, in particular, of electroweak processes in the 
Standard Model. These Sudakov-type corrections, which are found to exponentiate, 
are crucial for the high precision studies planned at the Next Linear (e"'"e~) Collider. 
Our approach is based on the use of an evolution equation in the infrared cut-off 
parameter fi, which in turn is based on a generalization of a gauge invariant dispersive 
method for photon Bremsstrahlung originally due to Gribov. 
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We have assumed that the energy of the process is much larger than all the masses 
in the theory. However New Physics could appear in the TeV region and it will modify 
our DL asymptotic forms. Because our approach is based on the infrared evolution 
equation, we need simply to change the initial conditions in this region to incorporate 
the new particles and/or new interactions of the existing particles. The infrared evolu- 
tion equation has the form of a renormalization group equation in the two-dimensional 
impact parameter subspace. 

The usual BFKL equation can be also considered as a simplified version of the 
renormalization group equation but in two-dimensional longitudinal subspace. In this 
case, in the leading logarithmic approximation ^f^logs ~ 1, instead of their dependence 
from the infrared cut-off fi, we have conformal invariance of the t-channel partial waves 
in impact parameter space. However, in the next-to-leading approximation, the effect 



of the running coupling leads to violation of conformal invariance [|T5|. After the break- 
ing of conformal invariance the equation takes the form of a quantized Callan-Symanzik 
equation [|16[. Therefore, even in the case of the BFKL equation, the renormalization 



group has its peculiarities. It is related to the fact, that the equation determines 
not only the anomalous dimensions of the operators of different twists, but also their 
relative contributions. 

The expressions for the DL asymptotics that we have obtained in the Standard 
Model can be used in perturbation theory to verify the first order, and to predict 
the higher order, expressions for scattering amplitudes, as has been done for the DL 
asymptotics in QCD |13|, |T^. Finally, we reiterate that the accurate calculation of 



scattering amplitudes is important because the effects of New Physics can be rather 
small. 



Added Note 

After this work was completed a paper by P. Ciafaloni and D. Comelli [|T^ appeared, in 
which the electroweak DL corrections are considered in the particular case of massless 
// production by a source which is a singlet with respect to the SU{2) x U{1) gauge 
group. The results of this paper (in particular, nonexponentiating DL corrections) are 
not in agreement with our corresponding results. The difference appears due to the 
fact that the method in [|18| consists of factorizing off the virtual gauge boson with 
the softest (i.e. with smallest frequency) momentum k^^ by computing external line 
insertions only, and in iterating this procedure by setting the virtual momentum k = 
in the left-over diagram. We believe that this approach, which is in disagreement with 
the method based on the infrared evolution equation, is not valid. 
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